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Abstract 

We discuss the lepton sector of a realistic string-inspired model based on the Pati- 
Salam S77(4) x SU(2)l X SU(2)r gauge group supplemented by a U(l) family sym- 
metry. The model involves third family Yukawa unification, predicts large tan/5 ~ 
50, and describes all fermion masses and mixing angles, including approximate bi- 
maximal mixing in the neutrino sector. Atmospheric neutrino mixing is achieved 
via a large 23 entry in the neutrino Yukawa matrix which can have important phe- 
nomenological effects. We find that the recent BNL result on the muon (g — 2) can be 
easily accommodated in a large portion of the SUSY parameter space of this model. 
Over this region of parameter space the model predicts a CP-even Higgs mass near 
115 GeV, and a rate for r — ► \i + 7 which is close to its current experimental limit. 
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Recently the BNL E821 Muon g-2 Collaboration has reported a precise measure- 
ment of the muon anomalous magnetic moment [|TJ a^exp) = (g — 2)/2, 

a^(exp) = (11, 659, 202 ± 15) x lO" 10 . (1) 

When combined with the other four most recent measurements the world average of 
a M is now higher than the Standard Model (SM) prediction, 

a^SM) = (11, 659, 160 ±7) x 1(T 10 (2) 

by (43 ± 16) x 10~ 10 which corresponds to a discrepancy of 2.6o\ It is well known 
that Supersymmetry (SUSY) gives an additional contribution to a^SM) which is 
dominated by the chargino exchange diagram and approximately given by 

where ai is the SU(2) gauge coupling, \i is the SUSY Higgs mass parameter, M 2 is 
SU (2) gaugino mass, is the muon mass, Msusy represents the heaviest sparticle 
mass in the loop, and tan/? is the ratio of Higgs vacuum expectation values (VEVs). 
Note that the sign of Aa^(SUSY) depends on the sign of \i (relative to M 2 ). 

Well before the experimental result from BNL was published, it was realised that 
the additional SUSY contribution Aa^(SUSY) could be of the correct order of magni- 
tude to be observed by E821 providing that tan f3 is sufficiently large, and the relevant 
superpartner masses Msusy are not too large @. Since the reported result, there 
has been a blizzard of theoretical papers, showing how the result may be accomo- 
dated within SUSY in detail and for various models 0. The general concensus of 
these recent studies is that numerically the additional SUSY contribution is sufficient 
to account for the discrepancy between the SM value and the experimental value, 
providing that tan j3 ~ 10 and Msusy ~ 500 GeV, and of course that the sign of /x is 
positive. 
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Large tan/? is also required in order to have a Higgs boson of mass 115 GeV 
||] (where the LEP signal has a significance of 2.9a) and it is encouraging that 
both signals point in the same direction of large tan (3. It is even more encour- 
aging that some well motivated unified models have long predicted that tan/3 is 
large. In particular models based on the gauge groups 50(10) or the Pati-Salam 
group 577(4) ® SU{2) L ® SU(2) R predict Yukawa unification which in turn implies 
tan/3 ~ 50 H Is experiment giving us a hint that Nature favours one of these 
Yukawa unification models which predict large tan (51 

There is a further piece of experimental evidence in favour of these models, namely 
that they both contain gauged SU(2)r symmetry and hence they both predict three 
right-handed neutrinos and hence non-zero neutrino masses. Thus in these models 
neutrino masses are compulsory, and not optional as in SU(5) for example. Su- 
perKamiokande evidence for atmospheric neutrino oscillations |7| has taught us that 
neutrino masses are non-zero and furthermore that the 23 mixing angle is almost 
maximal. The evidence for solar neutrino oscillations is almost as strong, although 
the conclusions are more ambiguous f| . A minimal interpretation of the atmospheric 
and solar data is to have a three neutrino hierarchy. A simple and natural interpreta- 
tion of the data is single right-handed neutrino dominance (SRHND) 0. In a large 
class of models, including those with SRHND, the large atmospheric mixing angle 
is due to large and equal couplings in the 23 and 33 entries of the Dirac neutrino 
Yukawa matrix (in the LR basis) 



Y v ~ 1 (4) 




corresponding to the dominant third right-handed neutrino coupling equally to the 
second and third lepton doublets. The see-saw mechanism yields a physical neutrino 
with a mass about 5 x 10 -2 eV consistent with the SuperKamiokande observation 
providing the third right-handed neutrino mass is Mr 3 ~ 3 x 10 14 GeV 13. 



The large off-diagonal Yukawa coupling in Eq.|] will have an important effect on 
the 23 block of the slepton doublet soft mass squared matrix m 2 L , when the renormal- 
isation group equations (RGEs) are run down from Mqut to the mass scale of the 
third right-handed neutrino M^ 3 . In order to see this it is instructive to examine the 
RGEs for m 2 L , 

dm\ I dm 2 L \ 
dt { dt ) Yv=Q 

- ^kY>£ + m^Yj + 2nm^^ (5) 
62ti z 1 J 

where m 2 N , m 2 Hu are the soft mass squareds of the right-handed sneutrinos and up-type 
Higgs doublet, A v is the soft trilinear mass parameter associated with the neutrino 
Yukawa coupling, and t = \b.(Mq UT / fi 2 ) , where /i is the MS scale. The first term on 
the right-hand side represents terms which do not depend on the neutrino Yukawa 
coupling. Assuming universal soft parameters at Mqut, m i(0) = m%(0) = m^I, 
where / is the unit matrix, and A u (0) = AY U , we have 

dm 2 L ( dm 2 L \ (3mg + A 2 



dt y dt J Y =0 16n 2 L u v 

where in the basis in which the charged lepton Yukawa couplings are diagonal, the 
first term on the right-hand side is diagonal. In running the RGEs between Mqut 
and M#3 the neutrino Yukawa couplings lead to an approximate contribution to the 
slepton mass squared matrix of 

5m 2 L - ~ ^ (^F) (Sml + A 2 ) [Y v Yj] * -0.1(3mg + A 2 ) [y u YJ] (7) 
Using the SRHND form of the neutrino Yukawa matrix in Eq.|] we find 

y„y! - 



Y,Y,] 



(6) 




and according to Eq.|^ the large neutrino Yukawa coupling in the 23 position will 
imply an off-diagonal 23 flavour violation in the slepton mass squared matrix which 



will be of order 5-10% of the diagonal soft mass squareds, and will be observable in 
the lepton flavour violating (LFV) process r — > fi + 7 . In addition the 22 entry of 
the slepton mass squared matrix will receive a 5-10% correction which again is due to 
the large 23 neutrino Yukawa coupling, and is much larger than the usual correction 
due to the diagonal muon Yukawa coupling which is very small. The large 22 entry 
in Eq|| will thus give a significant correction to the relation between the GUT scale 
soft mass parameters and the muon (g-2) estimates. The main purpose of the present 
paper is to explore these observable effects in the framework of a particular model 
which predicts Yukawa unification, and hence large tan /3, namely the string-inspired 



Pati-Salam model based on the gauge group S77(4) <g> SU(2) L <g> SU(2) R |L0 . 

For completeness we briefly review the string-inspired Pati-Salam model. As in 
50(10) the presence of the gauged SU(2) R predicts the existence of three right- 
handed neutrinos. However, unlike 50(10), there is no Higgs doublet-triplet splitting 
problem since both Higgs doublets are unified into a single multiplet h. Heavy Higgs 
H, H are introduced in order to break the symmetry. The model leads to third 
family Yukawa unification, as in minimal 5O(10), and the phenomenology of this was 
recently discusssed [|] . Although the Pati-Salam gauge group is not unified at the field 
theory level, it readily emerges from string constructions either in the perturbative 
fermionic constructions ||11]| , or in the more recent type I string constructions ||12|| , 
unlike 5O(10) which typically requires large Higgs representations which do not arise 
from the simplest string constructions. 

The Pati-Salam gauge group [K|, supplemented by a U(l) family symmetry, is 

SU(4)®SU(2) L ®SU(2) R ®U(1) (9) 

with left (L) and right (R) handed fermions transforming as Fi ~ (4, 2, 1) and F R ~ 
(4, 1,2) in the superfield multiplets 

^=fr^e-V do) 



L,R 



The Higgs h contains the two MSSM Higgs doublets and transforms as h ~ (1, 2, 2) 



(ii) 



The Higgs H,H transform as H ~ (4, 1,2), ~ (4, 1,2) and develop VEVs which 
break the Pati-Salam group, while 9, 9 are Pati-Salam singlets and develop VEVs 
which break the U(l) family symmetry. 

H,H=( J J J VH ),■■■ (12) 
^ d H d H d H e H J 

We assume for convenience that all symmetry breaking scales are at the GUT scale, 

< H >=< H >=< vh >~ M ~ 10 w GeV (13) 
< >=< 6 >~ M ~ 10 16 Ge1/ (14) 
The fermion mass operators (responsible for Yukawa matrices Y U) Y dl Y ei Yv) are 



^M^J(^J (i5) 

The third family is assumed to have zero U (1) charge, and the 33 operator is assumed 
to be the renormalisable operator with n = p = leading to Yukawa unification. 
The remaining operators have n > with varying group contractions involving HH 
leading to different Clebsch factors. The latter are responsible for vertical mass 
splittings within a generation. The mass splittings between different generations are 
described by operators with p > arising from different U(l) charge assignments to 



the different families. The Majorana mass operators (responsible for Mrr) are |T3 

i 



(PP)( HH ) ( HS \ ( d ) 



(16) 



We recently discussed [JT4] neutrino masses and mixing angles in the above string- 
inspired Pati-Salam model supplemented by a U{1) flavour symmetry. We used the 
SRHND mechanism, which may be implemented in the 422 model by having a 23 



operator with p = and n = 1 where the Clebsch is non-zero in the neutrino direction, 
but zero for charged fermions. This results in a natural explanation for atmospheric 
neutrinos via a hierarchical mass spectrum. We specifically focused on the LMA 
MSW solution since this is slightly preferred by the most recent fits, and assuming 
this a particular model of high energy Yukawa matrices which gave a good fit to 



all quark and lepton masses and mixing angles was discussed ||14|| . The numerical 
values of the high energy Yukawa matrices in this example are reproduced in Table I. 
To study lepton flavour violation focusing on the effects of the large off-diagonal 23 
entry in Y u , in this study we have further suppressed the tiny entries Y el2 , Y el3 , and 



Y ul3 compared to the values quoted in [14 . Note that with the suppression above 



the branching ratio BR([i — > e + 7 ) stays well below the experimental limit, without 
substantially changing the predictions of fermion masses and mixing angles. This 
demonstrates that this channel is more model dependent than r — > fx + 7 which is 
our main focus in this paper. 

The neutrino Yukawa matrix in Table I has a similar structure to that discussed 
in Eq.§ and has large approximately equal 23 and 33 elements. Thus the Yukawa 
matrices in Table I are examples of the effect that leads to 5-10% corrections to the 
23 block of the slepton mass squared matrix itll that we discussed previously. We 
now turn to a numerical discussion of these effects. 
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Table I. Yukawa matrices at Mqut (from ref.|Q) where the matrix elements of Mrr are in GeV. 

In our numerical analysis we have adopted a complete top-down approach JT5J . 
At the GUT scale we kept l/a GUT = 24.5223, M GUT = 3.0278 x 10 16 GeV, e 3 = 
(as(MouT) — oigut) I otcuT = —4.0568%, and the matrices in table I as fixed. Here 
otcuT = Oi2L = «i, and «3 = «4. For simplicity, the soft scalar masses of the 
MSSM superfields were introduced at the same scale. Including the D terms from 
the breaking of the Pati-Salam gauge group they read || 

m\ h + g\ D 2 

- {9l - 2^) D 2 



m, 



m 
m 

m 
m 
m 
rn 



m, 



ni 
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H,, 



m\ L -Zg\D 2 

m\ R + (Sgj - 2g 2 R ) D 2 

™ Fr + (3g 2 + 2g 2 R ) D 2 

m\-2g 2 R D 2 

m 2 h + 2g 2 R D 2 . 



(17) 



In the numerical analysis we kept the equality between the two soft SUSY breaking 



scalar masses mp L = mp R = mp. Two-loop RGEs for the dimensionless couplings 
and one-loop RGEs for the dimensionful couplings were used to run all couplings 
down to the scale M 3 # where the heaviest right-handed neutrino decoupled from the 
RGEs. Similar steps were taken for the lighter M 2 r and Mm scales, and finally with 
all three right-handed neutrinos decoupled the solutions for the MSSM couplings were 



computed at the Z scale, m^ and D in Eqs.|17| were varied to optimize radiative elec- 
troweak symmetry breaking (REWSB), which was checked at one loop following the 
effective potential method in [16| . As tan j3 determines the Higgs bilinear parameter 
Bfi, there is a redundancy in our procedure since two input parameters, rrih and D, 
determine one condition for the Higgs VEV of 246 GeV. This freedom was removed 
by favouring solutions with low CP odd Higgs mass m^o as a result of the observa- 
tion that values of m^o at the upper end of the range allowed by REWSB at a given 
(mp, M1/2) point are correlated, through the choice of the D, with low values for the 
stau mass which then in turn push the branching ratio BR(r — > /j, + 7 ) above the 
experimental limit. For this reason we introduced a mild penalty \ 2 into our analysis 
to favour REWSB solutions with low values of m^o. This top-down approach en- 
abled us to control the /i parameter as well as tan f3. We explored regions with \i low 
(/i = 120GeV) and high (/1 = 300GeV) Q. As a reference point we kept tan/3 = 50, 
and the universal trilinear coupling A = 0. An experimental lower bound on each 
sparticle mass was imposed. In particular, the most constraining are: the LEP limits 
on the charged SUSY masses (m^±,rrif > 105GeV), the CDF limit on the mass of 
the CP odd Higgs state (m A o > 105-110 GeV, valid for tan/? « 50) 0, and the 

requirement that the lightest SUSY particle should be neutral. 

1 For tan/3 as large as 50, fi ^> 300GeV leads to too large SUSY threshold corrections to the 
masses of the third generation fermions r and b unless the sparticles in the loop have masses well 
above the 1 TeV region. @ || 

2 Note that in this study we are primarily concerned with the lepton sector of the model and the 
effects of the large 23 element of the Y v in Eq.4. For this reason we drop two important constraints 
in the quark sector from the analysis. In particular, we do not consider the constraints imposed by 
the BR(b — > S7) and accept the b quark mass heavier than the value in [|l9| by about 15%. We 
assume that the complete theory at the high energy scale will induce additional corrections to the 
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The results are presented as plots in the (rrif, M1/2) plane. In figure [I] we show 
the best fit values for the quantities at the GUT scale which were varied to obtain the 
electroweak symmetry breaking. As explained in the previous paragraph these values 
are not unique, but preferred. We note that, clearly, the D terms in Eqs.|T7| are just a 
fraction of the scalar mass mp while the scalar higgs mass parameter mh is generally 
found to be greater than mp. The sharp turns in the contour lines of constant D 
below M1/2 ~ 4:00GeV result from the pseudoscalar Higgs mass m^o reaching the 
experimental lower bound, as demonstrated on plots (c) and (d) in figure 2. The 
parameters and D can still adapt to this change for Mi/ 2 < AOOGeV. The allowed 
(mj?, M1/2) region is finally bounded from below because of the too low chargino 
mass. This bound is at M 1/2 « 280GeV for /1 = 120GeV, and M 1/2 « 140GeV for 
\x = 300GeV. The region to the left of the contour lines is disallowed due to the stau 
lighter than any neutral SUSY particle. 

In figure ^| we plot the spectrum of the two neutral Higgs bosons h° and A . For 
low M1/2 their masses are degenerate while for higher values of My 2 the pseudoscalar 
Higgs becomes degenerate with the heavier of the two CP even Higgs states. Our 
analysis shows that the mass of the lighter CP even state is preferred to be in the 
range 112-117GeV for soft SUSY masses below 1 TeV. The pseudoscalar mass is 
quite sensitive to the magnitude of the D terms and, as was explained earlier, it was 
mildly pushed towards lower values as an additional condition on top of the REWSB 
conditions. 

Figure ^| represents the main results of this study. It shows that the constraints 
from the recent BNL experiment are consistent with all other constraints imposed 
on the model. In fact, as shown in plots (a) and (b) the BNL 2a region practically 

overlaps with the portion of the (mp, My 2) plane below 1 TeV allowed by the direct 

quark yukawa couplings possibly through a set of higher dimensional operators of the form (15) 
modifying the quark input parameters in table I. 
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sparticle searches. As promised in the text after Eq.[8] we also focused on the contri- 
bution to from the 22 entry in the slepton matrix in (|7|) generated by the large 23 
entry in (^J). In our numerical analysis the y 2 minimization procedure was extended 
to maximise this contribution. Nevertheless the maximum enhancement we found 
was on the level of 6%. 

The large 23 entry in (Q) makes an important contribution to the lepton flavour 
violating decay r — > ji + 7 . Plots (c) and (d) present the contour lines obtained in 
the same analysis. The computed values should be compared to the experimental 
upper limit BR(t — > [A + 7 ) < 1.1 x 1CT 6 at 90% C.L.. These predictions are quite 
robust. In order to reduce this branching ratio below the experimental limit over the 
entire plane we found we had to vary all initial parameters to rather extreme values, 
including lowering tan/3 as much as by 10 and increasing the trilinear parameter A 
into the TeV range. 

In conclusion, we have discussed the lepton sector of a realistic string-inspired 
model based on the Pati-Salam SU (4) x SU(2)l x SU(2)r gauge group supplemented 
by a U(l) family symmetry. The model involves third family Yukawa unification, pre- 
dicts large tan f3 ~ 50, and describes all fermion masses and mixing angles, including 
approximate bi-maximal mixing in the neutrino sector. In particular atmospheric 
neutrino mixing is achieved via a large 23 entry in the neutrino Yukawa matrix which 
we have shown to have important phenomenological effects. We find that the recent 
BNL result on the muon (g — 2) can be easily accommodated in a large portion of the 
SUSY parameter space of the model. Over this region of parameter space the model 
predicts a CP-even Higgs boson mass near 115 GeV, and a rate for r — > // + 7 which 
is close to the current experimental limit. We find it encouraging that all of these 
phenomenological features can be simultaneously accomodated within a simple string- 
inspired model such as the one considered in this study. 
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Figure 1: Contour lines of GUT scale parameters rrih and D determined by the 
condition of the radiative electroweak symmetry breaking, for two different values of 
fi. (a) m h , for fi = 120GeV. (b) m h , for fi = 300GeV. (c) D, for fi = 120GeV. (d) 
D, for n = 300GeV. Values in the plots are in GeV. In all plots tan/5 = 50, A = 0, 
and mp = mp L — Tn f r . 
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Figure 2: Contour lines of the light CP even Higgs mass m^o and pseudoscalar 
Higgs mass m^o, for two different values of /i. (a) m h o, for fi = 120GeV. (b) m h o, 
for /i = 300GeV. (c) m A o, for \i = 120GeV. (d) m^o, for \i = 300GeV. Values in 
the plots are in GeV. In all plots tan (3 = 50, A = 0, and m F = mp L = rriF R . 
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Figure 3: Contour lines of Sa^SUSY) x 10 10 and BR(r -> (i + 7 ) x 10 6 , for 
two different values of fx. (a) da^SUSY) x 10 10 , for (i = 120GeV. (b) 
5a^SUSY) x 10 10 , for /j, = 300GeV. (c) BR(r -> /1 + 7 )xl0 6 , for p = 120GeV. 
(d) BR(t — > /i + 7 )xl0 6 , for = 300GeV. In (a) and (b) the long-dashed curve 
marks the central value for a M not accounted for by the Standard Model, while the 
short-dashed curves mark the 2a limits of this quantity. The experimental upper 
limit on BR(t -»• fj, + 7 ) is 1.1 x 10 6 0. In all plots tan/5 = 50, A = 0, and 
mp = mp L = mp R . 
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